Modified embedded-atom method interatomic potentials for the Mg-Al alloy system 
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We developed new modified embedded-atom method (MEAM) interatomic potentials for the Mg- 
Al alloy system using a first-principles method based on density functional theory (DFT). The 
materials parameters, such as the cohesive energy, equilibrium atomic volume, and bulk modulus, 
were used to determine the MEAM parameters. Face-centered cubic, hexagonal close packed, and 
cubic rock salt structures were used as the reference structures for Al, Mg, and MgAl, respectively. 
The applicability of the new MEAM potentials to atomistic simulations for investigating Mg-Al 
alloys was demonstrated by performing simulations on Mg and Al atoms in a variety of geometries. 
The new MEAM potentials were used to calculate the adsorption energies of Al and Mg atoms on 
Al (111) and Mg (0001) surfaces. The formation energies and geometries of various point defects, 
such as vacancies, interstitial defects and substitutional defects, were also calculated. We found that 
the new MEAM potentials give a better overall agreement with DFT calculations and experiments 
when compared against the previously published MEAM potentials. 

PACS numbers: 61.50.Lt, 62.20.Dc, 61.72.Ji, 68.35.-p 
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I. INTRODUCTION 

Magnesium alloys are becoming increasingly impor- 
tant in many technological areas, including aerospace 
and automotive industries. The usage of magnesium die 
castings, for example, is increasing in the automotive 
industryiiSii 3 - due to the lower mass densities of magne- 
sium alloys compared with steel and aluminum, higher 
temperature capabilities and improved crash-worthiness 
over plastics. The primary magnesium alloys for die- 
casting applications are the magnesium-aluminum alloys 
such as AM50 and AM60Bi££. 

To meet the industrial demand for high-strength light- 
weight magnesium alloys, it is essential to obtain detailed 
understanding of the effect of individual alloying elements 
on the properties of magnesium alloys, especially among 
the main constituent elements, Mg and Al. The alloying 
elements can form interstitial or substitutional defects, or 
can precipitate into small particles creating complex in- 
terface structures. The interactions between these alloy- 
ing elements need to be investigated using atomistic sim- 
ulation techniques such as molecular dynamics or Monte 
Carlo simulations. These atomistic simulations require 
accurate atomic interaction potentials to compute the 
total energy of the system. First-principles calculations 
certainly can provide the most reliable interatomic po- 
tentials. However, realistic simulations of alloy systems 
often require a number of atoms that renders these meth- 
ods impractical - they either require too much computer 
memory or take too long to be completed in a reason- 
able amount of time. One alternative is to use (semi- 
empirical interaction potentials that can be evaluated 
efficiently, so that the atomistic approaches that use them 
can, in certain cases, handle systems with more than a 
million atoms. 



There are two additional essential features that are 
expected from a useful semiempirical approach besides 
its efficiency: reliability and flexibility. A reliable in- 
teratomic potential would accurately reproduce various 
fundamental physical properties of the relevant element 
or alloy, such as elastic, structural, and thermal prop- 
erties. Reliability also includes transferability. A trans- 
ferable interatomic potential would perform reasonably 
well even under circumstances that were not used during 
its construction phase. A flexible semiempirical approach 
can represent interaction potentials among a wide variety 
of elements and their alloys using a common mathemat- 
ical formalism. The modified embedded-atom method 
(MEAM) potential proposed by Baskes et al. was the 
first semiempirical atomic potential using a single for- 
malism for fee, bec, hep, diamond-structured materi- 
als and even gaseous elements, in good agreement with 
experiments or first-principles calculationsZ^&i^. The 
MEAM is an extension of the embedded-atom method 
(EAM)^ to include angular forces. The EAM was 
able to reproduce physical properties of many metals and 
impurities. The EAM was applied to hydrogen embrit- 
tlement in nickel^, and to nickel and palladium with 
hydrogen 1 ^. Cherne et al. made a careful comparison of 
MEAM and EAM calculations in a liquid nickel system^ 4 -. 

Atomistic simulations of a wide range of elements and 
alloys have been performed using the MEAM potentials. 
Baskes^ first proposed the MEAM method to obtain re- 
alistic shear behavior for silicon. Baskes et al.— then 
provided the MEAM model of silicon, germanium and 
their alloys. The MEAM was also applied to 26 single 
elements^ and to silicon-nickel alloys and interfaces 1 ^. 
Gall et alM used the MEAM to model the tensile 
debonding of an aluminum-silicon interface. Lee and 
Baskes — improved the MEAM to account for the second 
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nearest-neighbor interactions. Also, Huang et al.— used 
the MEAM and two other potentials to determine de- 
fect energetics in beta-SiC. The MEAM parameters for a 
nickel and molybdenum-silicon system were determined 
by Baske&i^'— . Recently, an effort has been made by 
Lee et al.— to create the MEAM potentials for Cu, Ag, 
Au, Ni, Pd, Pt, Al, and Pb, based on the first and the 
second nearest-neighbor MEAM. A new analytic modi- 
fied embedded-atom method (AMEAM) many-body po- 
tential was also proposed and applied to 17 hep metals, 
including M g 22 i 23 . For the Mg-Al alloy system, a set of 
EAM potentials have been developed using the "force 
matching" method by Liu et al.— . The structural prop- 
erties of various polytypes of carbon were described using 
a MEAM potential 2 ^. Finally, Potirniche et al.— used 
the MEAM to analyze damage evolution in a single crys- 
tal nickel. 

The purpose of the present work is to develop the 
MEAM potentials for aluminum, magnesium, and their 
alloy systems based on first-principles calculations us- 
ing density- functional theory (DFT). Energy calcula- 
tions and geometry optimizations of various structures 
were performed within the local-density approximation 
(LDA) 2 !^ using ultrasoft pseudopotential a 29 ' 30 ^. The 
cross pair potential was constructed by fitting elastic 
properties from DFT calculations for aluminum and mag- 
nesium in the Bl reference structure. First, the equilib- 
rium lattice parameter, cohesive energy, bulk modulus, 
trigonal and tetragonal shear moduli were determined 
from DFT calculations. The pair potential was then con- 
structed to fit the equilibrium volume and bulk modu- 
lus from ab initio calculations. Moreover, an effort has 
been made to match the sign of trigonal and tetragonal 
shear moduli. The new MEAM potentials were used to 
find the most energetically favorable structures for sin- 
gle elements and their pair combinations. The resulting 
energy- volume curves reasonably match the ab initio cal- 
culations. Satisfactory agreement of vacancy formation 
and stacking fault energies from DFT and MEAM cal- 
culations was found. Throughout this paper, the perfor- 
mance of our new potentials will be compared with the 
previously published MEAM potential a 21 ' 22 i 24 . 

The paper is organized in the following manner. In 
Sec. HH we give a brief review of the MEAM. In Sec. IIIII 
the procedure for determination of the MEAM parame- 
ters is presented along with the new MEAM interatomic 
potential parameters. Validation of the newly developed 
MEAM potentials is presented in Sec. IIV1 Different bulk 
structures, surface defects, and point defects calculations 
were performed and compared with DFT calculations 
and experiments. Finally, in Sec.fVl we discuss and sum- 
marize the results. 



II. MEAM THEORY 

The total energy £ of a system of atoms in the 
MEAM 32 is approximated as the sum of the atomic en- 



ergies 

i 

The energy of atom i consists of the embedding energy 
and the pair potential terms: 

E i = F i (p i ) + ^J2^( r io)- ( 2 ) 

F is the embedding function, pi is the background elec- 
tron density at the site of atom i, and <pij (rij) is the pair 
potential between atoms i and j separated by a distance 
r^. The embedding energy Fi {pi) represents the energy 
cost to insert atom i at a site where the background elec- 
tron density is pi. The embedding energy is given in the 
form 

F i (p i ) = A i E?p i ln(p i ), (3) 

where the sublimation energy Ef and parameter Ai de- 
pend on the element type of atom i. The background 
electron density pi is given by 

(o) 

ft = V( r .). ( 4 ) 

Pi 

where 

3 / (fc)\ 2 

fe=l \p\ J 

and 

G(r) = VTTT. (6) 

The zeroth and higher order densities, pf*\ , p\ , and 
Pi are given in Eqs. ([9]). The composition-dependent 
electron density scaling p® is given by 

P ° - p i0 Z i0 G (Tf) , (7) 

where pio is an element-dependent density scaling, Zio 
is the first nearest-neighbor coordination of the reference 
system, and Tf 1 is given by 

rr^iEfM^ (8) 
z « fe=1 

(k) 

where s\ is the shape factor that depends on the ref- 
erence structure for atom i. Shape factors for various 
structures are specified in the work of Baskcs 9 . The par- 
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tial electron densities are given by 
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where ry a is the a component of the displacement vector 
from atom i to atom j. Sij is the screening function 
between atoms i and j and is defined in Eqs. (Tl6|) . The 
atomic electron densities are computed as 



a(k) i \ 



(10) 



The screening function Sij is designed so that Sij = 1 
if atoms i and j are unscreened and within the cutoff 
radius r c , and = if they are completely screened or 
outside the cutoff radius. It varies smoothly between 
and 1 for partial screening. The total screening function 
is the product of a radial cutoff function and three body 
terms involving all other atoms in the system: 



Sij Sij fc 



Sikj — fc 



Ar 

n Si^ 

Cikj Cmin,ikj 



c, 



+/ 



1 + 2 



- r ij r ik + r ij T jk ~ r 



[ r ik r2 jk) 



fc{x)={ [1-{1-X) 4 



X>1 

< x < 1 
x < 



(16a) 
(16b) 

(16c) 
(16d) 

(16e) 



Note that C m - ln and C max can be defined separately for 
each i-j-k triplet, based on their element types. The 
parameter Ar controls the distance over which the radial 
cutoff is smoothed from 1 to near r = r c . 



where r® is the nearest-neighbor distance in the 
single-element reference structure and 0^ is element- 
dependent parameter. Finally, the average weighting fac- 
tors are given by 
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where tg fe j is an element-dependent parameter. 
The pair potential is given by 



<Aij ( r ij ) — 4>ij ( r ij ) Sij 
1 
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E^(r i j) = -E i j (1 + 4 ( rij )) e-W 




(11) 



(12) 

(13) 

(14) 
(15) 



where Eij, ay and r£ are element-dependent parame- 
ters and Zij depends upon the structure of the reference 
system. The background densities pi (rij ) in Eq. (|13p are 
the densities for the reference structure computed with 
interatomic spacing ry. 



III. DETERMINATION OF THE MEAM 
POTENTIAL PARAMETERS 



A. MEAM potentials for pure Al and Mg system 

The previously published MEAM parameters for Al^i 
and Mg£ served as the basis for the present work. For 
some of the surface and point defect calculations, how- 
ever, we observed less than satisfactory agreement be- 
tween DFT and the MEAM calculations when these orig- 
inal MEAM potential parameters were used. We fol- 
lowed similar procedures prescribed by Lee et al.— and 
Baskes^ to fine-tune the parameters and improve the 
overall agreement with experiments and DFT calcula- 
tions. For each element, several material parameters ob- 
tained from the reference structure are utilized to deter- 
mine the model parameters. These materials parame- 
ters include the cohesive energy, equilibrium atomic vol- 
ume, bulk modulus, and several elastic constants. The 
most stable crystal structures were chosen as the refer- 
ence structures, namely a face-centered cubic (fee) struc- 
ture for Al and a hexagonal close packed (hep) structure 
for Mg. The new parameters obtained from the present 
work are listed in Table |TJ The merits of these new poten- 
tials are demonstrated in various calculations described 
in Sec. ITVl 
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TABLE I: Set of the MEAM potential parameters for pure Al and Mg. E c is the cohesive energy, ao is the equilibrium lattice 
parameter, A is the scaling factor for the embedding energy, a is the exponential decay factor for the universal energy, /3' 0-3 - 1 
are the exponential decay factors for the atomic densities, t'- 0-3 -' are the weighting factors for the atomic densities, C max and 
Cmin are the screening parameters. The reference structures for Al and Mg are fee and hep, respectively. 



Element 


E c [oV] 


o Q [A] 


A 


a 


/?«» 


13^ 


13^ 


(5^ 


t (.o) 




t m 


t (3) 


Cmax 




Al 


3.353 


4.05 


1.07 


4.64 


2.04 


1.50 


6.0 


1.50 


1.00 


4.00 


-2.30 


8.01 


2.8 


2.0 


Mg 


1.55 


3.20 


1.11 


5.45 


2.70 


0.0 


0.35 


3.0 


1.00 


8.00 


4.10 


-2.00 


2.8 


2.0 



TABLE II: The MEAM potential parameters for the Mg-Al 
alloy system. E c is the cohesive energy, r e is the equilibrium 
nearest neighbor distance, a is exponential decay factor for 

the universal 6IlGrgy, Cmax cind Cmin 

are screening parameters, 

po is the density scaling factor. 



0.4575 



0.6 



Parameter 


Value 


E c [eV] 


(£ C A1 + E c MS )/2 


r«[A] 


2.821 


a 


4.915 


C min (Al-Mg-Al) 


0.0 


C min (Mg-Al-Mg) 


2.0 


C min (Al-Al-Mg) 


2.0 


C* mi n(Al-Mg-Mg) 


2.0 


C max (Al-Mg-Al) 


2.8 


C max (Mg-Al-Mg) 


2.8 


C* max (Al-Al-Mg) 


2.8 


C max (Al-Mg-Mg) 


2.8 


Po(Al) 


1.0 


Po(Mg) 


0.6 



> 




fee (DFT) + 

•' fee (MEAM) 

* hep (DFT) a 

hep (MEAM) 

bee (DFT) * 

bee (MEAM) 

sc (DFT) □ 

sc (MEAM) 

i I i L 

1.4 1. 



FIG. 1: Atomic energies (total energies per atom) as a func- 
tion of the atomic volume (volume per atom) for Al atoms 
in fee, hep, bee and simple cubic (sc) crystal structures. The 
energies are measured from the equilibrium atomic energy of 
fee structure. Volumes are scaled by the equilibrium atomic 
volume of the fee structure V fcc . 



B. MEAM potential for the Mg-Al alloy system 



The parameters of the MEAM potential for the Mg-Al 
alloy system were determined from a procedure similar to 
the one prescribed by Lee^. The parameters were con- 
structed to fit the elastic properties obtained from the 
DFT calculations for MgAl in the rock-salt (Bl) struc- 
ture, which was chosen to be the reference structure. The 
new parameters obtained from the present work are listed 
in Table [TTJ Primary emphasis was put on matching the 
equilibrium volume and the bulk modulus, which were 
reproduced exactly (see Table . 



IV. VALIDATION OF MEAM POTENTIALS 



We demonstrate the validity and the transferability of 
the new MEAM potentials by performing simulations on 
Al and Mg atoms in a variety of structural arrangements. 



A. Bulk 

1. Pure Al and Mg system 

To test the validity of the MEAM potentials for sin- 
gle elements, each element was put into fee, hep, body- 
centered cubic (bec), and simple cubic (sc) crystal struc- 
tures. The atomic energies for several atomic volumes 
near equilibrium atomic volume were calculated. The 
results were compared with those of DFT calculations, 
as shown in Fig. [TJ As expected, the curve for the fee 
structure produced by the MEAM potential retraces the 
results of DFT calculations nearly perfectly since fee was 
used as the reference structure during the potential con- 
struction process. The agreement between the MEAM 
potential and DFT for the hep structure is also remark- 
able. The most important result, however, is the fact that 
the new MEAM potential correctly identified fee as the 
most stable structure for Al. Furthermore, the sequence 
of the structures is correctly predicted in the order of sta- 
bility by the new Al MEAM potential. The relative co- 
hesive energies, with respect to the one for the fee struc- 
ture, are also in good agreement with the DFT calcula- 
tions, although the result for the simple cubic structure is 
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hep (DFT) 
hep (MEAM) 
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FIG. 2: Atomic energies of Mg as a function of the atomic 
volume in fee, hep, and bee cubic crystal structures. The 
energies are measured from the equilibrium atomic energy of 
the hep structure. Volumes are scaled by the equilibrium 
atomic volume of the hep structure V hcp . 



slightly underestimated. The relative equilibrium atomic 
volumes, with respect to the one for the fee structure, 
are also well reproduced. We point out that the equilib- 
rium atomic volume for fee Al, obtained by the MEAM 
potential (16.6lA 3 ), is slightly different from the one pre- 
dicted by DFT (15.76A 3 ). This is due to the fact that 
the MEAM parameters are fitted to reproduce the exper- 
imental volume, while DFT within LDA tends to under- 
estimate the equilibrium lattice constants by roughly 1% 
(see, e.g. Ref. I34I). 

Fig. [2] shows the atomic energy plot for Mg atoms in 
different crystal structures compared with the results of 
the DFT calculations. The hep structure was used as the 
reference structure for the Mg MEAM potential, and the 
DFT data points for this structure are accurately repro- 
duced. The sequence of the structures is again predicted 
correctly in the order of stability by the new Mg MEAM 
potential. The relative atomic energies, with respect to 
the one for the hep structure, are also in good agree- 
ment with the DFT calculations. Note that the scale of 
the vertical axis of Fig. [5] is six times larger than that of 
Fig. [TJ and that the largest error in relative atomic ener- 
gies (fee case) is in the order of 0.01 eV. Similar to the 
Al potential, equilibrium atomic volume for hep Mg in 
the MEAM is set to the experimental value of 23.16A 3 , 
while DFT predicts a smaller value of 21.54A 3 . Both the 
MEAM and the DFT methods prefer a c/a ratio close to 
0.994 of the ideal c/a ratio. 



2. Mg-Al alloy system 

To compare Mg-Al alloy systems with different stoi- 
chiometric coefficients, we define the heat of formation 



TABLE III: Elastic parameters for MgAl in a Bl structure 
from the MEAM and DFT calculations. The units of the 
heat of formation per atom Hi and the equilibrium atomic 
volume Vo are eV and A 3 , respectively. The units of the bulk 
modulus Bo and all elastic constants dj are GPa. 



Method H f 


V 


Bo C44 (Cn 


-C X2 )/2 


DFT 0.4575 
MEAM 0.4575 


22.4 
22.4 


38.4 -39.1 
38.4 -14.3 


40.6 
29.8 


per atom as 








H t = 


Etot — 




(17) 




N Mg + N A i 



where E to t is the total energy of the system, A^Mg and 
Na\ are the numbers of Mg and Al atoms in the system, 
£Mg and £ai are the total energies per atom for Mg and Al 
in their ideal bulk structures, respectively. Table HTT1 lists 
some of the material properties the new MEAM poten- 
tial reproduces for a MgAl compound in a Bl structure 
compared with the predictions of the DFT calculations. 
Due to our emphasis on the first three properties during 
the construction process, the last two columns show some 
discrepancies between the MEAM and DFT results. 

Fig. [3J shows the heat of formation per atom Hf for the 
Bl, B2 and B3 structures compared with the results from 
the DFT calculations. The Bl (cubic rock salt) struc- 
ture was used as the reference structure for the Mg-Al 
alloy MEAM potential. Fig. [3J shows that the reference 
structure is not the most stable structure in Mg-Al bi- 
nary systems. Again, the sequence of the structures in 
the order of stability is predicted correctly by the new 
MEAM potential for the Mg-Al alloy system. The rel- 
ative cohesive energies, with respect to the one for the 
reference structure, are also in good agreement with the 
DFT calculations. The equilibrium atomic volume and 
bulk modulus of Mg-Al in the Bl structure are repro- 
duced almost exactly. Note that the abscissa of the plot 
in Fig. [3] is the actual volume instead of the volume ratio 
used in Fig. Q] and Fig. 

To further demonstrate the validity of our new poten- 
tials, we also computed the heat of formation per atom 
for many intermetallic phases of Mg-Al alloys. The to- 
tal energy values in Eq. Q7| of Bl, B2, B3, CI, C3, C9, 
C15, DO3, DO9, A15, LI2 and A12 structures were eval- 
uated at the optimal atomic volume for each structure. 
The results from the MEAM calculations, compared with 
the ones from the DFT calculations, are summarized in 
Fig. 2J Although the Mg and Al atoms in these inter- 
metallic phases are in a chemical environment very dif- 
ferent from the one in the reference structure (Bl), the 
agreement between MEAM and DFT is quite satisfac- 
tory. In most cases, MEAM preserves the order of sta- 
bility predicted by DFT. The differences in the heat of 
formation per atom from MEAM and DFT are less than 
0.5 eV at most. However, we note that the MEAM failed 
to predict that the formation of one of the experimen- 
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FIG. 3: The heat of formation per atom for MgAl alloys in 
the Bl, B2 and B3 crystal structures. 



TABLE IV: Surface formation energies for fee Al and hep 
Mg. The units are mj/m 2 . The second column indicates if 
the structure was relaxed. Comparisons with other previously 
developed MEAM potentials are also given. 



Surface 


Relaxed 


MEAM" 


DFT 




Others 


Ref. 21 Ref. 22 Ref. 24° 


Al(lll) 


No 


737 


992 




913 


Al(lll) 


Yes 


731 


988 


629 


912 


Al(110) 


No 


1068 


1371 




1113 


Al(110) 


Yes 


1035 


1349 


948 


1107 


Al(100) 


No 


1025 


1213 




1012 


Al(100) 


Yes 


1025 


1212 


848 


1002 


Mg(0001) 


No 


604 


638 




500 


Mg(0001) 


Yes 


595 


637 




310 499 


Mg(1010) 


No 


642 


855 




517 


Mg(1010) 


Yes 


523 


846 




316 515 



"The present MEAM potential. 

'Calculated using EAM parameters extracted from Ref. \24 



> 

CD 




surface area E SU1 .{ is defined as 



40 50 60 
Atomic Percent Mg 



80 



FIG. 4: The heat of formation per atom for Mg-Al alloys in 
various intermetallic phases with different stoichiometric coef- 
ficients. The results obtained from the new Mg-Al MEAM po- 
tentials are compared with the DFT calculations. The struc- 
ture names are written next to the symbols (open triangles 
for DFT and crosses for MEAM). 



E smf = (E tot - Ne) /A, 



(18) 



where E to t is the total energy of the structure with a 
surface, N is the number of atoms in the structure, e 
is the total energy per atom in the bulk, and A is the 
surface area. Table IIVI shows the surface formation en- 
ergies of many different surfaces constructed from fee Al 
and hep Mg crystals. Results from the present MEAM 
potentials are in good agreement with the DFT calcu- 
lations, representing a significant improvement over two 
of the previously published MEAM potentials". How- 
ever, our MEAM potentials and the EAM potentials by 
Liu et al.— exhibit comparable levels of validity: our new 
MEAM potentials perform better for Mg surfaces while 
Liu's EAM potentials give better agreement for Al sur- 
faces. 



2. Stacking fault energies 



tally observed Mg-Al alloy structures", 7(Mgi 7 Ali 2 ) 
denoted as the A12 structure in Fig. [U as an exothermic 
process. In comparison, our DFT calculation correctly 
predicted the Hf for this structure to be a negative value 
(-0.017 eV). 



B. Surfaces 

1. Surface formation energies 

Semi-infinite surface is one of the simplest forms of 
defects. To test the transferability of the new MEAM 
potentials, surface formation energies for several different 
surfaces are computed. Surface formation energy per unit 



Stacking fault is another kind of structure that occurs 
frequently in real materials and provides a good test en- 
vironment for newly developed semiempirical potentials. 
Stacking fault energy per unit area is defined by 

E si = {E tot - Ne) /A, (19) 

where E to t is the total energy of the structure with a 
stacking fault, N is the number of atoms in the system, 
e is the total energy per atom in the bulk, and A is the 
unit cell area that is perpendicular to the stacking fault. 

For Al, three stacking fault types from Hirth and 
Lothe^i were examined and the results are listed in Ta- 
ble |Vl For the case of stacking fault type /, our MEAM 
result is in a good agreement with the available ex- 
perimental value, even though our DFT result is lower 
than the experimental value. In all cases considered, the 
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TABLE V: Stacking fault energies for Al. Results from the 
present MEAM and DFT calculations are compared. Stacking 
fault energies per unit area are given in mj/m 2 . 



Element 


Fault 


Relaxed 


MEAM (DFT) b 
Present Other" p ' 


Al 


I 


No 


150 (136) 


169 


Al 


I 


Yes 


146 (133) 


142 140-160 


Al 


E 


No 


150 (135) 


169 


Al 


E 


Yes 


148 (133) 


154 


Al 


T 


No 


75 (62) 


84 


Al 


T 


Yes 


74 (61) 


77 



"Calculated using EAM parameters extracted from Ref. |24 
'Experimental results from Ref. l37l. 

/ = ABCBCABC 

E = ABC ABC B ABC ABC 

T = ABCABCABACBACB 



TABLE VII: Adsorption energies E^as and optimized height 
above the surface from MEAM calculations. Results from 
DFT calculations are given in parentheses. Units are eV and 
A for energies and heights, respectively. 



Surface 


Adatom (site) 


-Eads 


Height 


Al(lll) 
Al(lll) 
Al(lll) 
Al(lll) 


Al (hep) 
Al (fee) 
Mg (hep) 
Mg (fee) 


-2.64 (-3.29) 
-2.67 (-3.26) 
-1.70 (-1.09) 
-1.70 (-1.07) 


2.19 (2.05) 
2.17 (2.08) 

2.46 (2.35) 

2.47 (2.35) 


Mg(0001) 
Mg(0001) 
Mg(0001) 
Mg(0001) 


Al (hep) 
Al (fee) 
Mg (hep) 
Mg (fee) 


-2.17 (-2.68) 
-2.17 (-2.68) 
-1.43 (-0.81) 
-1.49 (-0.82) 


2.11 (2.16) 
2.09 (2.14) 
2.17 (2.28) 
2.37 (2.27) 



Adsorption on surfaces 



TABLE VI: Stacking fault energies for Mg. Results from the 
present MEAM and DFT calculations are compared. Stacking 
fault energies per unit area are given in mj/m 2 . Comparisons 
with other previously developed MEAM potentials are also 
given. 



Element 


Fault 


MEAM" 


EAM 6 


AMEAM C 


DFT d 


DFT C 


Mg 


h 


7 


27 


4 


(18) 


(20.9) 


Mg 


h 


15 


54 


8 


(37) 


(43.7) 


Mg 


T 2 


15 


54 




(45) 


(51.3) 


Mg 


E 


22 


81 


12 


(61) 


(68.1) 



"The present MEAM potential. 

'Calculated using EAM parameters extracted from Ref. l24i 

C AMEAM results in Ref. H| 

"DFT results from the present study. 

E DFT results in Ref. HI 

h = ABABABCBCBCB 

I 2 = ABABABCACACB 

T 2 = ABABABCBABAB 

E = ABABABCABABAB 



present MEAM potential shows better overall agreements 
with DFT calculations compared with the EAM poten- 
tial by Liu et al.— . 

For Mg, four stacking fault types from the calcula- 
tion of Chetty and Weinert^ were examined. Total en- 
ergy calculations for ii, I2, T%, and E stacking fault 
types were performed using both DFT and MEAM cal- 
culations. The results are compared in Table IVIl The 
present MEAM potential shows a substantial improve- 
ment over the previously published MEAM potential by 
Hu et al.— . The stacking fault energies are consistently 
underestimated by the present MEAM potentials com- 
pared to the results of the DFT calculations, while the 
results by the EAM potential from Ref. [3 are consis- 
tently overestimated. 



The adsorption energy of a single adatom E a ds is given 

by 

E a ds = E tot — E SUI f — E atom , (20) 

where E tot is the total energy of the structure with the 
adatom adsorbed on the surface, E sur f is the total energy 
of the surface without the adatom, and -Eatom is the to- 
tal energy of an isolated atom. We placed single Al and 
Mg atoms at the fee and hep sites (see, e.g., L0vvik and 
Olsen— for the definition of these sites) on Al(lll) and 
Mg(0001) surfaces. The entire structures were then re- 
laxed to determine the adsorption energies. The height 
of the relaxed adatom is measured from the farthest atom 
(the one least affected by the adsorption) in the top sur- 
face layer. Table IViTl shows the results obtained from the 
MEAM and DFT calculations. Even though there are 
small quantitative discrepancies, the qualitative agree- 
ment between these two sets of data is quite satisfactory. 
For instance, DFT calculations predict that on Al(lll) 
surfaces the adsorption of the Al atoms (same kind) is 
stronger (bigger adsorption energies and shorter heights) 
than the adsorption of Mg atoms (different kind). On the 
other hand, DFT calculations predict that on Mg(0001) 
surfaces, the adsorption of the Mg atoms (same kind) is 
weaker (smaller adsorption energies and longer heights) 
than the adsorption of Al atoms (different kind). Both 
of these features are clearly demonstrated by the new 
MEAM potentials. 

C. Point defects 

1. Vacancy 

The formation energy of a single vacancy EJ ac is de- 
fined as the energy cost to create a vacancy: 

Er c = E tot [N]-Ne, (21) 

where E tot [N] is the total energy of a system with N 
atoms containing a vacancy and e is the energy per atom 
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TABLE VIII: Calculated single vacancy properties. Single va- 
cancy formation energy EJ^ and formation volume Q v values 
are obtained from the relaxed structures containing single va- 
cancies. Here Qo is the bulk atomic volume. All energy values 
are listed in eV. The results from the ME AM calculations are 
compared with the results from the DFT calculations given 
inside the parentheses. 



Present Others Present Others 

Al 0.68 (0.67) 0.68°, 0.68" 0.66 (0.76) 0.72°, 0.61" 
Mg 0.58 (0.82) 0.59 c , 0.87° 0.76 (0.75) 0.83 c , 0.88" 

"MEAM results in Ref. 

'Calculated using EAM parameters extracted from Ref. l24i 
C AMEAM results in Ref. 



TABLE IX: The formation energies of various kinds of inter- 
stitial point defects in Al and Mg. All energy values are given 
in eV. The results from the MEAM calculations are compared 
with the results from the DFT calculations given inside the 
parentheses. 



Bulk (structure) 


Interstitial (site) 


MEAM (DFT) 


Al (fee) 


Al (dumbbell) 


2.32 (2.94) 


Al (fee) 


Al (octahedral) 


2.91 (3.06) 


Al (fee) 


Al (tetrahedral) 


3.14 (3.68) 


Al (fee) 


Mg (octahedral) 


2.77 (3.79) 


Al (fee) 


Mg (tetrahedral) 


5.09 (4.25) 


Mg (hep) 


Mg (octahedral) 


1.29 (2.36) 


Mg (hep) 


Mg (tetrahedral) 


1.53 (2.35) 


Mg (hep) 


Al (octahedral) 


2.13 (1.97) 


Mg (hep) 


Al (tetrahedral) 


2.79 (2.11) 



in the bulk. Table I Villi shows the formation energy of 
single vacancies for fee Al and hep Mg obtained from the 
MEAM and DFT calculations. The new MEAM poten- 
tials reproduced a DFT value for Al vacancy formation 
energy very well, although the value for Mg was esti- 
mated somewhat low. Furthermore, the present MEAM 
potentials reproduce the correct amount of reduction in 
volume due to the formation of a vacancy. This also 
represents a substantial improvement over the existing 
MEAM potentials. 



2. Interstitial point defects 

The formation energy of an interstitial point defect 
E l f nt is given by 

Eg* = E tot [N + A] - E tot [N] - e A (22) 

where E tQt [N] is the total energy of a system with N (Mg 
or Al) atoms, E to t [N + A] is the total energy of a system 
with N atoms plus one atom of type-^4 (Mg or Al) in- 
serted at one of the interstitial sites, and ea is the total 
energy per atom of type->l in its most stable bulk struc- 
ture. Note that the inserted atom A can be the same type 



TABLE X: The formation energies of substitutional point de- 
fects in Al and Mg. All energy values are given in eV. 



Bulk (structure) 


Substitute 


MEAM (DFT) 


Al (fee) 


Mg 


0.38 (0.06) 


Mg (hep) 


Al 


0.55 (0.09) 



as the matrix, in which case the point defect becomes a 
so-called self- interstitial defect. Interstitial atom forma- 
tion energies were calculated for Al and Mg at octahedral, 
tetrahedral, and dumbbell sites. Atomic position and 
volume relaxation were performed. The results of these 
calculations arc listed in Table HXl to be compared with 
the results from the DFT calculations. DFT results are 
well reproduced in general. According to the present cal- 
culations, the most stable form of a self- interstitial defect 
for fee Al crystal is a dumbbell along the [100] direction, 
in agreement with the DFT results and an experimental 
observation by Jesson et al.— . The new MEAM poten- 
tials, however, failed to reproduce the results of the DFT 
calculations for the self-interstitial defects in a hep Mg 
crystal. Our new MEAM potential indicates that the 
octahedral site will be more stable than the tetrahedral 
site, while the DFT calculations predict that both sites 
will have nearly the same formation energies. In both 
of the heterogeneous interstitial defects, the new MEAM 
potentials produce the same relative stability of different 
interstitial sites with the DFT calculations. However, 
Table IIXI should not be used to predict the most stable 
interstitial defects; it was not the purpose of the present 
work to perform an exhaustive search to draw such con- 
clusions. 



3. Substitutional point defects 

The formation energy of a substitutional point defect 
Ef uh , in the case of the substitution of an Al atom with 
a Mg atom, is defined by 

E! uh = E tot [Mg A1 ] - E tot [AIaj] - e Mg + 6ai (23) 

where £t t[Mg A1 ] is the total energy of a system of Al 
atoms plus one Mg atom that replaced an Al atom, 
i?tot[AlAi] is the total energy of the original system of 
Al atoms without a defect, £M g and £ai are the total en- 
ergies per atom for Mg and Al in their ideal bulk struc- 
tures. The formation energy of a substitutional point 
defect for other cases can be defined similarly. Table IXl 
shows the results of substitutional defect calculations us- 
ing the MEAM potentials and the DFT method. The 
new MEAM potentials predict correctly that substitut- 
ing a Mg atom in a hep structure with an Al atom costs 
more energy than the reverse as indicated by the DFT 
results, although the formation energies in both cases are 
larger than the values from the DFT calculations. 
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D. Molecular dynamics simulations 

To validate the new potentials for molecular dynam- 
ics simulations, we calculated the melting temperatures 
of pure Al and Mg crystals. We followed the procedure 
prescribed by Morris et al.— to establish co-existence of 
solid and liquid phases to determine the melting temper- 
atures. We obtained 930 K for pure Al crystal, which 
is in excellent agreement with the experimental value of 
933 K. For Mg crystal, however, the two-phase method 
did not give a satisfactory result: 530 K compared to the 
experimental value of 923 K. To compare with other po- 
tentials, we followed a single-phase method as described 
by Kim and Tomanek— , in which the temperature is in- 
creased at a constant rate and the specific heat of the 
system is monitored. Using this method, we obtained 
780 K as the melting temperature of Mg crystals. This 
result is comparable to 745 K obtained by Liu et al.— 
using an EAM potential and following a similar method 
to compute the melting temperature. 

The difficulty in computing accurate melting temper- 
atures for hep metals, such as Mg, using semiempirical 
potentials is well documented^. We believe that this is 
related to a small energy difference between hep and fee 
structures in our MEAM potential for Mg (see Fig. \2§ 
and instability of hep structures in MEAM^. Further 
detailed investigation of this subject is beyond the scope 
of the present work and will be reported in separate pa- 
pers. 

V. CONCLUSIONS 

In this study we developed a new set of MEAM po- 
tentials for Al, Mg and their alloy systems using first- 



principles calculations based on DFT. The validity and 
transferability of the new MEAM potentials were tested 
rigorously by calculating physical properties of the Mg- 
Al alloy systems in many different atomic arrangements 
such as bulk, surface, and point defect structures. The 
new MEAM potentials show a significant improvement 
over the previously published potentials, especially for 
the surface formation, stacking faults, and point defect 
calculations. The new Mg-Al alloy potentials, however, 
failed to predict the stability of the 7(Mgi7Alx2) alloy 
intermetallic phase, suggesting the need for further im- 
provements. 
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